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O ■ ABSTRACT 

o 

O We demonstrate the importance of near-infrared radiation from hot dust for Comp- 

ton cooling of electrons/positrons in quasar jets. In our model, we assume that the 
non-thermal radiation spectra observed in OVV quasars are produced by relativistic 
electrons/positrons accelerated in thin shells which propagate down the jet with rel- 
ativistic speeds. We show that the Comptonization of the near-IR flux is likely to 
dominate the radiative output of OVV quasars in the energy range from tens of keV 



up to hundreds of MeV, where it exceeds that produced by Comptonization of the UV 
radiation reprocessed and rescattered in the Broad Emission Line (BEL) region. The 
main reason for this lies in the fact that the jet encounters the ambient IR radiation 
over a relatively large distance as compared to the distance where the energy density of 
the broad emission line light peaks. In the soft - to mid energy X-ray band, the spectral 
component resulting from Comptonization of the near-IR radiation joins smoothly with 
the synchrotron-self-Compton component, which may be responsible for the soft X-ray 
flux. At the highest observed 7-ray energies, in the GeV range, Comptonization of 
broad emission lines dominates over other components. 
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1. INTRODUCTION 

Recent discovery that blazars are strong and variable 7-ray emitters (von Montigny et al. 
1995) provided independent evidence that the blazar phenomenon (smooth continuum emission in 
all observable bands; large-amplitude, rapid variability; and high linear polarization) is produced 
by relativistic sub-parsec jets which are oriented at small angles to the line of sight (Dondi & 
Ghisellini 1995; for a recent review, see, e.g., Ulrich, Maraschi, & Urry 1997). Radio loud quasars 
with such an orientation of the jet are observed as optically-violent-variable (OVV) and highly 
polarized quasars, while the line-weak radio galaxies are seen as BL Lac objects. 

Non-thermal spectra of blazars can be divided into two components, the low energy com- 
ponent, commonly interpreted in terms of the synchrotron radiation mechanism, and the high 
energy component, presumably produced by the inverse-Compton (IC) process (see, e.g., review by 
Sikora 1997). In BL Lac objects, the IC energy losses of relativistic electrons/positrons are likely to 
be dominated by the synchrotron self-Compton (SSC) process (Mastichiadis & Kirk 1997; Ghisellini 
et al. 1998; Coppi & Aharonian 1999), whereas in the OVV quasars, Comptonization of the diffuse 
ambient radiation field is probably more important (Sikora, Begelman, &; Rees 1994 [hereafter: 
SBR94]; Blandford & Levinson 1995). The above, however, should be regarded as a trend rather 
than the rule, and, e.g., the case of BL Lacertae shows that Comptonization of external radiation 
can be important also in BL Lac objects (Madejski et al. 1999). 

The models proposing that 7-rays in OVV quasars result from Comptonization of the diffuse 
ambient radiation field are commonly known as external radiation Compton (ERC) models. There 
are many variants of the ERC models. In some of them, radiation sources are approximated by 
homogeneous "blobs" propagating along the jet (SBR94; Ghisellini et al. 1998; Kusunose, Takahara, 
& Li 2000), while in others — by an inhomogeneous flow (Blandford & Levinson 1995). They also 
differ regarding the dominant diffuse ambient radiation field. For small distances, < 10 16 ~ 17 cm, 
this can be provided directly by the accretion disc (Dermer & Schlickeiser 1993), while at larger 
distances, the broad emission lines (BEL) and near-IR radiation from hot dust are likely to dominate 
(SBR94). In the "blob" model, the distance along the jet where most of the 7-ray flux is produced 
can be estimated from variability time scale and from the location of the spectral break where the 
luminosity of the Compton component peaks. The 1 -j- 3 day time scales of 7-ray outbursts and 
typical location of their spectral breaks in the 1 4- 30 MeV range are consistent with production 
of 7-rays by Comptonization of light from BEL regions and dust at distances 10 17 — 10 18 cm (see 
Table 1 in SBR94). 

The situation in the X-ray band is more complex. In part, the SSC process may be important 
in the soft/mid energy X-ray bands (Inoue & Takahara 1996; Kubo et al. 1998). Furthermore, X- 
ray variability usually shows lower amplitude than variability in the 7-ray band (Sambruna 1997; 
Wehrle et al. 1998). This suggests that X-rays produced co-spatially with 7-rays are contami- 
nated by another source of X-rays, possibly located at larger distances along the jet and therefore 
varying on much longer time scales. Superposition of radiation components produced at lower and 
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higher distances can also explain progressively weaker variability in the synchrotron component at 
lower frequencies (Brown et al. 1989; Edelson 1992; Hartman et al. 1996; Ulrich, Maraschi & Urry 
1997). Since both sources located at smaller and larger distances are related to dissipation events 
presumably triggered by collisions between inhomogeneities propagating down the jet at different 
velocities, one can expect that location and strength of these collisions change with time. As a 
result, the relative amplitude of different spectral components may also be changing from epoch to 
epoch. The best studied example is the quasar 3C 279, which is usually very variable in the optical 
band (this is an OVV quasar), but during the high energy flare in 1996, the optical flux varied less 
than 20% (Wehrle et al. 1998). 

As was demonstrated by Sikora & Madejski (2000), studies of X-ray radiation in OVV quasars 
allow a verification of the presence of low energy electrons in a jet. This, in turn, can be used to 
estimate the pair content of quasar jets. We conduct such studies in this paper by more detailed 
modeling of production of radiation. This is also the first work where the contribution of all three 
Compton components: SSC, C(BEL) and C(IR) is taken into account simultaneously. 

This paper is organized as follows: in §2 we describe our model of non-thermal flare production 
in relativistic jets. In §3 we use the model to reproduce typical spectra observed in 7-ray detected 
OVV quasars and illustrate their dependence on the low energy break in the electron/positron 
injection function. We use two classes of models, one where the radiative output is dominated 
by Comptonization of IR radiation from hot dust, and another where 7-ray flux results from 
Comptonization of broad emission lines. In that section, we also discuss the model parameters and, 
in particular, calculate the pair content of the jet plasma, required to provide the observed spectra. 
We summarize our main conclusions in §4. 



In our model, we assume that relativistic electrons/positrons are enclosed within a very thin 
shell which propagates along the jet with relativistic speed. With this, we can describe evolution of 
the electron energy distribution, 7V 7 , using the following continuity equation (see, e.g., Moderski, 
Sikora, & Bulik 2000) 



2. 



THE MODEL 



2.1. Evolution of the Electron Energy Distribution 




(1) 



where the rate of electron/positron energy losses is 




3 r 



(2) 



r is the distance from the central engine, T is the bulk Lorentz factor, (5 = \/T 2 — 1/T, 7 is 
the random Lorentz factor of an electron/positron, Q is the rate of injection of relativistic elec- 
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trons/positrons, and dr = j3cTdt' . The second term on the right-hand side of Eq. (2) represents the 
adiabatic losses for two-dimensional expansion of the jet. Radiative energy losses are dominated 
by: 



(1) synchrotron radiation, 

u B l , (3) 



dt 1 J s 3m e c 
(2) Comptonization of synchrotron radiation, 

d<y\ 4<j T _ , _ o 



and (3) Comptonization of external radiation, 

U ) erc = ~ 9^ (U * //(flBL) + "*//(^) r T ' ^ 

where = B' 2 /8tt is the magnetic field energy density, u' s is the energy density of the synchrotron 
radiation field, ^//(bbl) — Lb ^Ji e lnr is the energy density of the broad emission line field present 
at the actual distance of a source/shell propagating downstream a jet, ^//(/r) — £lr ^<tsb T a /c 
is the energy density of the infrared radiation field, is the fraction of the central radiation 
reprocessed into near infrared by hot dust, and T is temperature of dust assumed to be located at 
r m > r (see Eq. 17). Note that the expression for Udiff(BEL) does not include the contribution 
from emission lines produced at much larger or smaller distances than r. This is because for any 
smooth radial distribution of external diffuse radiation field, as is likely to be the case of BELR in 
quasars (Peterson 1993), the lines produced over larger distances (> r) provide much lower energy 
densities, while the lines produced at smaller distances are very significantly redshifted in the shell 
frame. Here, as elsewhere in the paper, primed quantities denote measurements made in the source 
co-moving frame. Energy density of the synchrotron radiation is given by: 



U ' s ~ 2W^ 



f'^S.max 

/ ^L'g^dW (6) 

Jv' h 

abs 



where a is the cross section radius of the jet at a distance where the outburst is produced, u' abs is 
the frequency at which the optical thickness due to synchrotron-self- absorption is equal to unity 3 , 
hv' S ,max = (mi 2 ma AB'/B cr )m e c\ B cr ~ 4.4 x 10 13 Gauss, and 



d7 



dt' 



(7) 



It should be mentioned here that in the case of thin shell geometry, the synchrotron energy density, 
given by formula (6), is reached throughout the source on much shorter time scale than a/c, and, 



3 We calculate u' abs using an approximate formula given by SBR94 and derived assuming spherical geometry of the 
source. For the shell geometry, u' abs is larger and angle dependent, but never exceeds the value for spherical case by 
more than < (a/A') 2/7 . 
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therefore, unlike the case of spherical sources, the SSC flares are not expected to lag significantly 
behind the synchrotron flares. We also note that our assumption about the comoving thinness of 
shells is well justified by the detailed studies of shocks produced by collisions of inhomogeneities 
moving with different velocities, provided their Lorentz factors differ by less than a factor 2 (see, 
e.g., Komissarov and Falle 1997). 



2.2. Modeling of the Radiation Spectra 



Since we consider only the contribution to the overall spectra by the regions of the jet which 
are optically thin, the observed spectra can be computed using the following formula 



where 



"M*) - ft 11 v ' L 'v t r ; Q'P 4 d cos e d<p , 

c(3t 



(8) 



(9) 



1-/0COS0' 

V = [r(l — /3cos O)]^ 1 = T(l + (3 cos 6') is the Doppler factor, 6 is the angle between velocity of the 
shell element and direction to the observer, v = XV, and t is the observed time (see Appendix A). 
The intrinsic synchrotron luminosity, u'L' Sl>l , is given by Eq. (7). The intrinsic SSC luminosity is 
(see, e.g., Chiang & Dermer 1999) 



v ' L ' SSC a 



V3&T ,3/2 

8n jr 2 v 



I3u> 



T I 1-3/2 , / 

L sy v s dz/c 



where 



and 



Max 



3i/ 



^ I max 



Min 



v s,- 



3i/ 

A. ,2 

mm 



I mm J 



(10) 



(11) 



(12) 



In contrast to the synchrotron and SSC emission which is isotropic in the co-moving frame, 
the radiation produced by Comptonization of external radiation is anisotropic, and thus 



div'L',.. ,) i 



d7 



C(i) 



where 



d7 



dt' 



C(i) 



[9'} 



-7 V m, 



3m e c 



and i: BEL or IR. The right-hand side of Eq. (13) is calculated for 



(13) 



(14) 



(15) 
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Anisotropy of the ERC process (Eq. 14) and its astrophysical implications have been extensively 
discussed by Dermer (1995), and some general comments on the significance of this anisotropy can 
be found in Sikora (1997). 

2.3. Application to Homogeneous ERC Models 

Using the formalism given above, we consider homogeneous ERC models. Specifically, our 
goal is to compare two scenarios: in the first case, 7-rays are produced by Comptonization of 
near-IR radiation (which we subsequently denote as models A), while in the second case, it is due 
to Comptonization of broad emission line photons (hereafter models B). Our model assumptions 
are: 

• electron/positron injection function is a power-law, Q = Kj~ p , for 7& < 7 < imax, and has a low 

energy tail, Q oc 7 -1 , for 7 < 7^; 

• 76 is lower than the energy of the break, 7 C , produced due to cooling effect (SBR94 4 ) and, 

therefore, it is the latter which determines the position of luminosity peaks of the spectral 
components C(BEL) and C(IR); 

• electrons/positrons are injected at a constant rate within a distance range (ro; 2ro) and uniformly 

fill the shell whose radial width is A' < vq/T; 

• the shell propagates down the conical jet with a constant Lorentz factor T. The half-opening 

angle of the jet is Oj = 1/T; 

• the intensity of the magnetic field is B(r) = (ro/r)Bo, the energy density of broad emission lines 

is u<iiff(BEL) = ( r o/ r ) 2 ^B£;L,0) and the energy density of infrared radiation is uir = const; 

• the observer is located at an angle 6 b s = 1/T. 

With the above assumptions, the model approximates a situation where the shell contain- 
ing relativistic plasma is formed due to collision of two perturbances moving down the jet with 
different velocities (see, e.g., SBR94 and Rees & Meszaros 1994). The low energy break in the 
electron/positron injection function, 75, corresponds to the characteristic energy of pre-heated elec- 
trons/positrons (see, e.g., Kirk, Rieger, & Mastichiadis 1998), and the low energy tail below this 
break is introduced to mimic a limited efficiency of the pre-heating process, nature of which is very 
uncertain (Hoshino et al. 1992; Levinson 1996; McClements et al. 1997). 

Input parameters of the models were chosen to reproduce typical features of OVV quasar flares: 
location of the 7-ray luminosity peak in the 1-30 MeV range; apparent 7-ray luminosity in the 
range 10 47-48 erg s" 1 ; apparent synchrotron luminosity in the range 10 46 ~ 47 erg s" 1 ; deficiency of 
radiation in the soft/mid X-ray bands; riicixiiniim synchrotron rrGCjusncy, rnax ^ 

10 15 Hz; and 

time scale of 7-ray flares on the order of days. For broadband spectra of blazars see, e.g., von 
Montigny et al. (1995) and Fossati et al. (1998); for time scales of rapid (daily) 7-ray variability 



4 Note that in SBR94 the cooling energy break is denoted by 75, while in our paper, it denotes the low energy 
break in the injection function. 
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see, e.g., Michelson et al. (1994), Mattox et al. (1997) and Wehrle et al. (1998). 



3. RESULTS AND DISCUSSION 

As it was shown by SBR94, the high energy 7-rays detected by EGRET on board CGRO 
from OVV quasars can be produced by Comptonization of broad emission line light as well as by 
Comptonization of infrared radiation (see Table 1 in SBR94). Whereas the broad emission line 
flux gained a lot of attention as a good candidate for the dominant source of seed photons for 
the inverse Compton process in a jet, the infrared radiation from dust was largely ignored in the 
literature (see, however, Wagner et al. 1995a and references therein). To the best of our knowledge, 
this paper is the first where production of both spectral components, C(IR) and C(BEL), is treated 
simultaneously. We discuss the energy density of those components below. 

Regarding the diffuse radiation field on sub-parsec scales provided by broad emission line 
clouds, the line luminosities in radio loud quasars are typically in the range 10 44 — 10 46 erg s _1 
(Celotti, Padovani & Ghisellini 1997; Cao & Jiang 1999). Assuming that the scale of the broad 
emission line region is on the order of the distance at which non-thermal flares are produced in a 
jet, one can find that the corresponding energy density at a distance r = 10 18 rig cm is 

Lbel 1fl _3 Lbel,45 -1 /,^x 
u diff(BEL) - j^r c - 3 x 10 erg s , (16) 

where Lbel,45 = £REL/10 45 ergs _1 . For the diffuse IR radiation, just as in the cases of Seyfert 
galaxies and radio quiet quasars, the spectra of lobe-dominated radio loud quasars show very 
prominent near-IR bumps (Sanders et al. 1989; de Vries et al. 1998). These bumps are commonly 
interpreted as thermal radiation produced by hot dust. Such dust is expected to be heated by UV 
radiation of an accretion disc and located in the innermost parts of a geometrically thick molecular 
torus, at a distance 

/ L UV , ! n i8 L UV,46 /•-,--> 

r ^vi^^" 4x10 ^r cm ' (17) 

where Ljjy = 10 46 Luv,m erg s _1 is the luminosity of the accretion disc, and T = 10 3 T3 K is 
the effective temperature of dust. The dust provides radiation field with energy density which at 
f < Tir is approximately constant and equal to 

uir ~ Or ±°SB T 4 /c ~ 2.3 x 1(T 3 fe R /0.3) T 4 erg s" 1 . (18) 

The model parameters are given in Table 2, and illustrated in Figures 1-3, which show 
time-averaged broad-band spectra computed using the model equations and assumptions given in 
§2. Specifically, we consider three cases: (1) a pure C(IR) model (A, shown in Fig. 1); (2) pure 
C(BEL) model (B, shown in Fig. 2); and (3) a combination of the two (A+B, shown in Fig. 3). 
In all three figures, a comparison of the four panels illustrates the dependence of the time averaged 
spectra on the low energy break in the injection function, 75. 

To validate the models in greater detail, we also consider spectra obtained for the two well- 
studied objects. In Fig. 4 we present our model fits to the Jan-Feb 1996 outburst in 3C 279 (Wehrle 
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et al. 1998) and to the March 1995 outburst in PKS 0528+134 (Blazejowski et al. 1997; Kubo et al. 
1998), using a combination of models A and B. The model parameters for those fits are presented 
in Table 2. 



3.1. Spectra 

Figures 1-3 clearly illustrate that the general features of spectra of 7-ray detected OVV 
quasars - in particular, the high ratio of Compton to synchrotron luminosities, and the comparative 
deficit of X-ray radiation - can be reproduced by both types of models, A (where C(IR) dominates) 
and B (where C(BEL) dominates). However, the details of the spectra are different. The models 
A yield X-ray spectra which are much smoother and closer to the observed ones than the spectra 
produced by models B. This results from the fact that in order to reproduce the observed spectra, 
in models A, the Lorentz factors of electrons/positrons 7 C (where the cooling break occurs) need 
to be larger than in models B, causing the X-ray portion of the SSC component to be harder in 
models A as compared to B. 

Another important difference between these two models is in the high energy 7-ray spectral 
band. The spectral components due to C(IR) have a break at ~ 1 GeV, whereas EGRET/ CGR O 
observations show that blazar spectra extend at least up to 5 GeV. Therefore, the EGRET spectra 
can be reproduced only as a superposition of C(IR) and C(BEL) (see Fig. 3). This also provides 
an attractive explanation for 7-ray spectra being steeper than the X-ray spectra by more than the 
value predicted by cooling effect, i.e., Aa = 0.5. Examples of spectra of objects with Aa > 0.5, 
often called "MeV blazars," can be found in McNaron-Brown et al. (1995), Blom et al. (1995), and 
Malizia et al. (1999). 



3.2. Low Energy Injection Break and Pair Content 

Total number of relativistic electrons plus positrons injected in the shell during the flare is 

fjmax m K ( 1 \ 

N, = At'l Q , 7 = -£_ I (,„ 7t+ _), ( 19 ) 

where At' = Ar/cT = tq/cT is the injection time interval. The number of protons enclosed in the 
shell is on the order of 

N p ~ N p At x ~ , (20) 

m p c z T c 

where L p is the energy flux of cold protons in a jet and A is the width of the shell as measured in 
the external frame. We note here that because the energy dissipated in a jet in our model results 
from collisions between inhomogenities, one should not expect heating of protons to much larger 
average energies than mildly relativistic. Furthermore, this energy is quickly lost due to adiabatic 
expansion and/or converted back to the bulk energy. This is why jet energy flux is expected to be 
dominated by cold protons. 
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Combining Eqs. (19) and (20), we find 

N + N e m p c 2 r K 



In 7b + 



1 



) 



(21) 



N p 2N p 2L p XjP 



>-i 



p-1 



(We denned N e as a sum of numbers of both relativistic electrons and positrons, and thus for N + > 
a few xN p the charge conservation gives N e ~ 2N + .) Assuming the kinetic luminosity of the jet to 
be Lj = 10 47 Lj t 4j ~ 10 47 erg s _1 (see, e.g., Rawlings Sz Saunders 1991; Celotti & Fabian 1993) and 
noting that A < tq/T 2 (see assumption about the shell width in §2), we obtain results presented in 



The pair content implied by our fits to the time averaged spectrum of the outburst in 3C 279 
is n + /n p ~ 4, and for the outburst in PKS 0528+134 — n + /n p ~ 17. Note, however, that our fits 
are not unique, since such parameters as time scales of flares, energy densities of external radiation 
fields, and jet powers are not very well constrained by observations. 

As we can see from Table 3, the pair content ranges from a few for 7& ~ 100 up to more 
than hundred for 7b ~ 1. The recent detection of circular polarization in several compact radio 
sources (Wardle et al. 1998; Homan & Wardle 1999) and its interpretation in terms of the Faraday 
conversion process favor low values of 7b, and therefore a large pair content. However, there are 
some objects with such hard X-ray spectra (Malizia et al. 1999; Fabian et al. 1998), that they can 
only be explained if 7b 3> 1. But even in these cases, as long as 7b < 100, at least some pair content 
is required. Such pairs can be produced via interaction of a cold proto-jet with radiation of the hot 
accretion disc corona (Sikora & Madejski 2000). 



As we can see from Table 1, in the models A and A+B, uir ~ 10 _4 erg s _1 , i.e. about 20 times 
less than energy density of radiation produced by silicate grains with its maximum temperature 
~ 1000 K, and about 100 times less than energy density of radiation produced by graphite grains 
with its maximum temperature ~ 1500 K. This implies the value of uir to be lower than expected, 
in order to avoid overproduction of X-rays. In either model, the corresponding dust temperature 
is ~ 500 K, i.e. 2-3 times lower than maximum, and the distance of dust from the central source 
is ~ 5y/Ljjv,46 P c - Alternatively, overproduction of X-rays by Comptonization of IR radiation 
can be avoided by postulating that the 7-ray luminosity peak is not defined by the value of the 
cooling break j c as assumed above, but rather, by the injection break 7b (see Ghisellini et al. 1999; 
Mukherjee et al. 1999). 

In the model A+B, Udiff(BEL) — u diff(iR) — 10~ 4 er g cm~ 3 (see Table 1), which is 30 times 
lower than the average in radio loud quasars. However, as recent reverberation campaigns show, 
broad emission lines have peak luminosities at distances ~ 3 x 10 1 y/Luv,46 cm ( see review by 
Peterson 1993), and therefore the amount of the diffuse line radiation at distances where flares are 
produced (~ ro-^2ro ~ 10 18 cm) is likely to be a small fraction of the total Lbel- In case of models 
B, Udiff(BEL) is much closer to the observed values and there is no need for significant reduction 
of Lbel at distances corresponding to tq. 



Table 3. 



3.3. 



External Radiation Fields 
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3.4. Bulk Compton Radiation 

The models presented above don't take into account the so-called bulk-Compton process - 
the Comptonization of ambient diffuse radiation by cold electrons/positrons in a jet (Begelman & 
Sikora 1996; Sikora et al. 1997). This process should lead to the production of narrow bumps at 
energies T 2 hv ex t, i.e., ~ 2 x (r/15) 2 keV and ~ 0.3 x (r/15) 2 keV for v ext = vjjv an d v ext = vir, 
respectively. The largest contribution should come from very small distances, where the jet just 
starts to be relativistic and collimated. The absence of soft X-ray bumps in the observed spectra 
of OVV quasars (see, e.g., Lawson & M c Hardy 1998) suggests that the region of jet formation is 
very extended, possibly reaching distances not very much smaller than those where most of blazar 
radiation is produced, or else, the number of pairs in a jet is very small (cf. Sikora et al. 1997). 
Noting that the number flux of cold pairs at r < r$ should at least be equal to the number flux of 
relativistic pairs accelerated at r ~ vq, the latter condition corresponds to 7& 3> 1. 

3.5. Variability 

Fig. 5 shows the time evolution of the broad-band spectrum presented in Fig. 3, binned in 
the time intervals At = to = ro/2r 2 c. As one can see from those Figures, the flares decay at 
different rates in different spectral bands. Despite the transverse size of the source and related 
light travel effects, the dependence of energy losses on electron/positron energy is manifested very 
strongly, causing the synchrotron and Compton spectra above the luminosity peak to steepen with 
time. Additional steepening of the 7-ray spectrum is caused by the fact that C(BEL) drops faster 
with distance than C(IR). This is because the energy density of broad emission line light decreases 
with distance, while energy density of IR radiation is roughly constant. Such a steepening of 7- 
ray spectra is observationally confirmed both by direct observations of the slope changes during 
individual flares (Miicke et al. 1996; Mukherjee et al. 1996; Mattox et al. 1997), and statistically, 
by comparison of time averaged spectra as measured during different epochs (Pohl et al. 1997). 

Steepening of the synchrotron spectrum is also consistent with observations, which show that 
the amplitude of variability increases with the increasing energy of the observing band. However, 
in many OVV quasars, the amplitude of the variability of the synchrotron component, even at the 
highest frequencies, is much smaller than amplitude of variability in 7-rays. A particularly clear 
example of such behavior is provided by the simultaneous multi-band observations of the Jan-Feb 
1996 flare in 3C 279 (Wehrle et al. 1998). This can be explained by assuming that the synchrotron 
component is heavily contaminated by radiation produced at larger distances and/or by stationary 
shocks. 

The situation at the keV energies is quite complex, but very interesting, as it provides important 
constraints on models. If the spectra measured in this range are entirely dominated by ERC 
components, then the drop of flux is predicted to be achromatic (with constant slope). This 
is because the low energy tails of ERC components are produced by electrons /positrons whose 
energy losses are dominated by adiabatic expansion. However, as suggested by Inoue & Takahara 
(1996) and Kubo et al. (1998) and confirmed by our models, the low energy X-rays are likely to 
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be dominated by the SSC component. Then the whole X-ray spectrum is a superposition of a 
softer SSC component and harder ERC component, which should result in a concave shape. Such 
spectra are in fact in a few cases observed directly, and are also implied by the fact that two-point 
slopes between the X-ray and 7-ray bands are harder than the X-ray spectra measured in the 
soft/mid X-ray bands alone (Comastri et al. 1997). As we can see from Fig. 5, such X-ray spectra 
- formed by a superposition of SSC and ERC components - harden as the source fades. This is 
caused by the SSC component, with luminosity which is proportional to the square of synchrotron 
luminosity, decaying faster than the low energy tail of ERC components. Hardening of X-ray 
spectra during the decay of the flare was in fact observed in some objects (see Ghisellini et al. 1999; 
Malizia et al. 1999). In most cases, however, the X-ray spectrum softens as the flare fades (Lawson, 
M c Hardy & Marscher 1999; Comastri et al. 1997). This discrepancy can be resolved by postulating 
that the X-ray spectra, particularly at the lowest energies, are contaminated by radiation produced 
either at larger distances (see, e.g., Unwin et al. 1997), or in stationary shocks which may form if 
jet is reconfined and/or bent (see Komissarov & Falle 1997 and references therein). Furthermore, 
the light curves of outbursts can be further modified by a change of shell velocity (see, e.g., Dermer 
&; Chiang 1998), or by a change of direction of propagation. The latter case seems to be favored 
by VLBI observations, which often show strong curvature of jets in OVV quasars on parsec scales 
(Appl, Sol & Vicente 1996; Bower et al. 1997; Rantakyro et al. 1998; Tingay, Murphy & Edwards 
1998). Finally we emphasize that even very high amplitude flares are never isolated: observations 
show that they smoothly join with lower amplitude neighboring flares, forming with them longer 
lasting high states. In order to explain the continuously flaring high state light curves in terms 
of propagating shells of relativistic plasma, one must assume that number of shells per a distance 
range tq — 2tq is on the order of N s h ~ T 2 . For this particular number, the number of shells 
contributing to the observed radiation at a given instant is N s h, bs ~ N s h/T 2 ~ 1 (see Appendix in 
Sikora et al. 1997). For a much lower number of shells, the observer would see very isolated flares, 
while for much larger number of shells the resultant fluctuation (amplitude of the smaller flares 
superposed on the "high state") would be smaller than is commonly observed. 

4. CONCLUSIONS 

• Comptonization of infrared radiation produced by hot dust should be taken into account in all 

radiation models of blazars. Such radiation, likely to be present in quasar cores on parsec 
scales, is sufficiently dense to compete with broad emission line flux as the source of Compton 
cooling of relativistic electrons/positrons in parsec/sub-parsec jets; 

• Both types of ERC models, either those where the 7-ray luminosity is dominated by C(IR) 

(models A) or those where the 7-ray luminosity is dominated by C(BEL) (models B), are able 
to reproduce all basic broadband spectral features of OVV quasars during their outbursts; 

• If the blazar radiation is produced at distances larger than the distance at which the energy 

density of the broad emission line light peaks, the C(IR) component is expected to dominate 
over the C(BEL) component; 
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• The ERC models predict the general steepening of the 7-ray and synchrotron spectra with energy, 

in a qualitative agreement with observations. However, the lower observed variability amplitude 
of the synchrotron flux as compared with the 7-ray flux suggests a strong contamination of 
synchrotron radiation produced in the process of a flare by synchrotron radiation produced at 
larger distances or in stationary shocks; 

• X-ray spectra in OVV quasars consist of SSC and ERC components, and the latter can be either 

C(IR) or C(BEL). Just as in the case of synchrotron radiation, additional contribution to the 
soft/mid energy X-rays is expected from quasi-steady component produced by distant and/or 
stationary shocks. This can explain why the amplitude of variability in the X-ray bands is 
smaller than in 7-rays; 

• The ERC models of 7-ray production predict the pair content of the jet plasma, n + /n p , to be in 

the range from a few x £7,47 for 75 ~ 100 up to hundred x Lj^ for 7^ ~ 1. Low values of 75 
(and therefore a large e + /e~ pair content) are favored by observations of circular polarization, 
assuming its interpretation in terms of the Faraday conversion process is correct. 

The project was partially supported by Polish KBN grant 2P03D 00415, ITP/NSF grant 
PHY94-07194, NASA grants NAG-6337 and NASA observing/ADP grants to the University of 
Maryland and USRA. 
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A. RELATIVISTIC TRANSFORMATION OF OBSERVED LUMINOSITY 

(1 ) The case of a point source of radiation 

Observed monochromatic luminosity of a point source at the moment t is 

rlP HP' 

5LM = 4n l£ [t] = ^ = V * 5L 'ArAs\ , (Al) 

(Rybicki & Lightman 1979), where 

cBt , . , 

r = T^^. < A2 > 

v -W^r) sr{l+l}C0S<, ' h (A3) 

u' = f/T>, and 8 is the angle between the direction of the source motion and the direction to the 
observer, n b s - 

(2) The case of an extended source of radiation 

Since we are only interested in the optically thin radiation, the observed luminosity of an 
extended source can be computed by summing contribution from its small pieces, treated as point 
sources. Of course, because of light travel effects, for pieces moving along different 0's, contribution 
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must be taken from different distances (see Eq. A2). Using approximation of infinitely thin shell 
(A <C ro/r 2 ), we divide the source into 50, j = ilj/N pieces, where N must be enough large to allow 
treat individual pieces as point sources (i.e., with negligible light travel effects). The observed 
luminosity is then given by 

{vL v )[t\ = Sat {uSL u )[9;t] = Y, N V\v'8L' v ,)[r-n' obs ] . (A4) 

For uniform shells, one can introduce the formal quantity, "intrinsic luminosity" L' u [9; t] = SL' u [0; t] x 
N = 6L' v [9;t] x (Clj/SQj), and with this and iV — > oo, Eq. (A4) can be converted into Eq. (8). 
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FIGURE CAPTIONS 

Fig. 1. — Time averaged model spectra of blazars for four values of 75. The dotted line marks the 
typical slope of the X-ray spectra observed in OVV-type blazars (a = 0.6); the dashed lines enclose 
the 1-30 keV X-ray band. In each panel we show three spectral components: synchrotron (SYN), 
synchrotron-self-Compton (SSC) and Comptonization of near-IR dust radiation [C(IR)]. (Model 
A). (For the model parameters, see Table 1). 

Fig. 2.— Same as Fig. 1, but C(IR) is replaced with C(BEL) (Model B). 

Fig. 3.— Same as Fig. 1, where both C(IR) and C(BEL) are included (Model A+B). 

Fig. 4. — Left panel: model fit to the outburst spectrum observed in PKS 0528+134 (Model A+B). 
The X-ray spectra (ASCA) and 7-rays (EGRET) are simultaneous (McNaron-Brown et al. 1995; 
Kubo et al. 1998). Radio and optical data are nonsimultaneous and are taken from the archive 
sources (see Blazejowski et al. 1997 for details and references therein). Right panel: model fit to 
the outburst spectrum observed in 3C 279. Radio, optical, X-ray and 7-ray data are simultaneous 
(Wehrle et al. 1998) (Model A+B). 

Fig. 5. — Binned spectra in the time intervals At = ro/2r 2 c (Model A+B). 
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Table 1. Parameters used in models 



MODEL 


r 


ro [cm] 


B [G] 






lb 


'"/max 


P 






A 


15 


7.0 x 10 17 


0.43 


6.8 x 10 49 


1 


- 100 


6.0 x 10 3 


2.2 


0.0 


10- 4 


B 


15 


6.0 x 10 17 


1.34 


6.9 x 10 49 


1 


- 100 


1.0 x 10 4 


2.2 


2.3 x 10- 3 


0.0 


A+B 


15 


7.0 x 10 17 


0.43 


8.9 x 10 49 


1 


- 100 


6.0 x 10 3 


2.2 


10- 4 


10- 4 



Table 2. Parameters used in modeling 3C 279 and PKS 0528+134 



OBJECT 


r 


ro [cm] 


B [G] 




lb 


"fmax 


P 






3C 279 


20 


7.0 X 10 17 


0.81 


2.8 x 10 so 


150 


6.5 x 10 3 


2.4 


0.6 x 10 -3 


0.2 x 10 -3 


PKS 0528+134 


15 


7.0 x 10 17 


1.7 


1.2 x 10 so 


25 


4.7 x 10 3 


2.2 


2.1 x 10~ 3 


4.7 x 10- 4 



Table 3. {n + /n p ) m i n ratios. 



7b 


1 


5 


20 


100 


A 


96 


40 


12 


2 


B 


97 


41 


12 


2 
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FIG.2. 
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